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ABSTRACT 

We present simple, analytic solutions to the Einstein-Maxwell equation, which 
describe an arbitrary number of charged black holes in a spacetime with positive 
cosmological constant A. In the limit A = 0, these solutions reduce to the well 
known Majumdar-Papapetrou (MP) solutions. Like the MP solutions, each black 
hole in a A > solution has charge Q equal to its mass M, up to a possible overall 
sign. Unlike the A = limit, however, solutions with A > are highly dynamical. 
The black holes move with respect to one another, following natural trajectories 
in the background deSitter spacetime. Black holes moving apart eventually go 
out of causal contact. Black holes on approaching trajectories ultimately merge. 
To our knowledge, these solutions give the first analytic description of coalescing 
black holes. Likewise, the thermodynamics of the A > solutions is quite interest- 
ing. Taken individually, a \Q\ — M black hole is in thermal equilibrium with the 
background deSitter Hawking radiation. With more than one black hole, because 
the solutions are not static, no global equilibrium temperature can be defined. In 
appropriate limits, however, when the black holes are either close together or far 
apart, approximate equilibrium states are established. 



* This work was supported in part by NSF grant NSF-THY-8714-684-A01 



Introduction 



In this paper we give simple, analytic solutions to the Einstein-Maxwell equa- 
tion, which describe collections of charged black holes in a spacetime with positive 
cosmological constant A. In the limit of vanishing cosmological constant, our so- 
lutions reduce to the well known Majumdar-Papapetrou (MP) solutions [1]. For 
all values of A > 0, individual black holes in the solutions have charge Q equal 
to mass M, up to a possible overall sign. For A > 0, however, this is no longer 
the condition for extremality [2]. We study the mechanical and thermodynamic 
properties of the A > solutions, which turn out to generalize the properties of 
the MP solutions in interesting ways. 

For A = 0, the black holes are static. In contrast, the black holes in a A > 
solution are highly dynamical. The black holes ignore one another and follow nat- 
ural trajectories in the background deSitter spacetime. The black holes eventually 
either merge, or move out of causal contact. As far as we know, these solutions 
give the first analytic description of coalescing black holes. However, there are 
many questions about the causal structure and dynamics of these solutions which 
remain to be answered. 

The thermodynamics of solutions with A > is also quite interesting. The 
black holes in the MP solution each have vanishing Hawking temperature, and 
so are in thermal equilibrium, with each other and also with the background flat 
spacetime. With A > 0, the background deSitter spacetime has an ambient Hawk- 
ing temperature [3]. A single \Q\ — M black hole actually has a temperature equal 
to the deSitter temperature, and hence would be in a state of thermal equilibrium 
with the background spacetime. The temperature, however, depends on both the 
cosmological constant and the mass of the black hole. Adding more black holes, 
with differing masses, then does not produce a thermal equilibrium state. Indeed, 
with more than one black hole, because the solutions are not static, no global equi- 
librium temperature can be defined. Approximate temperatures, however, can be 
defined in the limits where the black holes are either widely separated or coalesced. 
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A given solution, then, describes a transition between one state of thermal equi- 
librium in the far past and another equilibrium state in the far future. However, 
an understanding of the nonequilibrium thermodynamics of the transition between 
these states will require further investigation. 

Majumdar-Papapetrou Solutions 

In Newtonian mechanics, a collection of charged point particles, each having 
charge equal to its mass , can stay at rest in a state of mechanical equilibrium. 
The electrostatic repulsions of the particles exactly balance the gravitational at- 
tractions. Remarkably, the same balance holds in general relativity. The MP 
solutions to the source free Einstein-Maxwell equation correspond to this Newto- 
nian situation. The MP solutions, themselves, are geodesically incomplete. Hartle 
and Hawking [4] showed how the MP solutions could be analytically extended and 
interpreted as a system of charged black holes. The metric and gauge field for the 
MP solutions are given by 

ds 2 = -n- 2 dt 2 + VI 2 (dx 2 + dy 2 + dz 2 ) , A t = ft -1 , 
Q=l + J2—, r i = J(x-x i )* + (y-y i )* + (z-z i )*, & 

~ n v 

where m\ and (xi, m,Zi) are the masses and locations of the black holes. It can be 
shown [4] that the points r-{ = actually represent event horizons of area 4irm 2 . For 
the case of one black hole, the metric (1) is just the extreme Reissner- Nordstrom 
metric in isotropic coordinates. 

A single Reissner-Nordstrom black hole having charge equal to its mass is the 
simplest example of an extremal black hole. If the charge of the hole were in- 
creased (or the mass decreased) further, then the curvature singularity would no 
longer be hidden behind an event horizon. If such a naked singularity arose in 

* We use geometrical units, G = c = h = 1 
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a 'physical' solution, it would violate the cosmic censorship conjectures and lead 
to a breakdown of predictability. That the Hawking temperature of an extremal 
Reissner-Nordstrom black hole vanishes is impressive evidence for cosmic censor- 
ship. The evaporation of a charged black hole terminates when it reaches the 
extremal state, leaving the singularity hidden^. The MP solutions then describe 
charged black holes in thermal, as well as mechanical, equilibrium with tempera- 
ture equal to zero. A related property is that the MP solutions are exact ground 
states of N = 2 supergravity. Since all quantum corrections to the effective action 
expanded around the MP solutions vanish [6] . Analogues of the MP solutions have 
also been written down for dilaton black holes [7,8]. The individual black holes in 
these solutions are also extremal ones. 



Reissner-Nordstrom-deSitter Solutions 

There is an analogue of the Reissner-Nordstrom solution for a spacetime with a 
cosmological constant. The Reissner-Nordstrom-deSitter (RNdS) metric and gauge 
field in Schwarzschild coordinates are given by 

ds 2 = -V{R)dT 2 + V~ 1 (R)dR 2 + R 2 dQ 2 , A T = -5 

(1) 

2M Q 2 1 , K J 

where M and Q are the mass and charge of the hole and A is the cosmological 
constant. The metric has a curvature singularity at the origin. There are event 
horizons at the radii where V(R) vanishes. Taking M — Q — in (1) gives the 
static form of the deSitter metric. 

Romans [2] has recently studied the thermodynamics of these solutions. This 
is more complicated, because for A > the deSitter Horizon also radiates at its 

f Holzhey and Wilczek [5] have studied how the evaporation process for charged dilaton black 
holes can terminate at an extremal state, even though the temperature of this state may 
be nonzero 
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own temperature [3]. Indeed, Gibbons and Hawking [3] have extended the laws of 
black hole thermodynamics to include cosmological event horizons. The Hawking 
temperature for a horizon at radius p is given by 




V\P)\, 



(2) 



where k is the surface gravity at the horizon. As in the A = case, Romans finds 
extreme RNdS black holes with zero temperature, in which the inner and outer 
black hole horizons coincide. For A > 0, the extremal black holes have M < \Q\. 
In Appendix A, we show that, as in the A = case [9], it is impossible to destroy 
the horizon of an RNdS black hole by throwing in charged particles to charge it 
past the extremal limit. 

Another interesting class of RNdS solutions have the temperatures of the black 
hole and deSitter horizons equal. Remarkably, these solutions have \Q\ — M*. The 
metric function then has the form 



In the naive picture of black hole evaporation, the \Q\ — M solutions are ther- 
modynamically stable states and are the end points of the evaporation process. If 
M > \Q\ then the black hole is hotter than the deSitter horizon and it will evap- 
orate until it reaches M — \Q\. If M < \Q\, then the deSitter horizon is hotter 
and the black hole will accrete radiation until it reaches M — \Q\. It is interesting 
that for A > the conditions of extremality and thermal equilibrium no longer 
coincide. 

* This had been noted previously by Mellor and Moss [10] 



V(R) = (1 - - ^ 



(3) 



The common temperature is given by [2] 
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We note that the total gravitational entropy of an RNdS black hole is maxi- 
mized, for fixed charge and cosmological constant, in this equilibrium state. The 
total gravitational entropy is given by the sum of the areas of the black hole and 
cosmological event horizons, 



The result then follows from the generalized first law of thermodynamics given in 
[3]. For an infinitesimal perturbation between RNdS solutions of fixed charge, the 
generalized first law states that 



where k^s and kbh are the surface gravities of the deSitter and black hole event 
horizons. The change in entropy of such a variation can then be written as 



Clearly, the entropy is extremized for k^s = —k-bh, which coincides with the 
condition that the black hole and deSitter temperatures be equal, and happens 
when \Q\ = M. From the discussion in the previous paragraph, we see that the 
extremum is a maximum. This result is in contrast to an evaporating Schwarschild 
black hole, where the gravitational entropy decreases. Here, the deSitter horizon 
acts like the walls of a box. The entropy of the deSitter horizon measures, in some 
sense, the entropy of the Hawking radiation carried away from, or absorbed by the 
black hole. 



Sgrav = -7 A BH + T^dS- 



(5) 



Kds8A d s - k-bhSAbh = 0, 



(6) 
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Motion of Test Particles 



Are cosmological analogues of the MP solutions built out of extremal or \Q\ = 
M black holes? The black holes in the MP solution "ignore" one another. We 
would like to find some similar phenomenon for charged black holes in a deSitter 
background. To find the right criterion, we look at the motion of charged test 
particles in the RNdS metric (1). We will see that the motion of q = m test 
particles in a Q = M RNdS background is particularly simple. 

The conserved energy £ of a test particle of charge q, mass m and 4-velocity 
u a on a radial path in an RNdS background is 

- = -e («« + ^~A a ) 

m \ m / /i\ 

dT a ^ ' 

= V(R)±--±A T (R), 
dr m 

where £ a is the static Killing vector and r is the proper time. Together with the 
normalization condition u a u a — — 1, this gives the equation of motion 



(-£) = -V'(fi) + + ^M«)) 2 . (2) 



\ dr / mm 

Substituting in the gauge field and metric function for a Q = M black hole, this 
becomes 

fdR\ 2 f, M\ 2 1,„, f, qM 



If the test particle has q = m and has energy to equal its rest mass (i.e. E = m), 
then this further reduces to 

dR /A„ 

— = ±\-R. 4) 
dr V 3 V ; 

This in turn is the same as the equation of motion for a minimum energy test 

particle in a background deSitter spacetime (i.e. Q — M — 0). This looks like 

what we want. The q = m test particle is, in some sense, not affected by the 

presence of the Q = M black hole. This hint will turn out to be what we need to 

guess an exact solution. 
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Note that there are two choices for the path of the q = m test particle in (4). 
It can be either "ingoing" or "outgoing". This is in contrast to the A = case, 
where a minimum energy q = m test particle stays at rest in the field of a Q = M 
black hole. Choosing one path or another breaks the time reversal invariance of 
the system. We will continue to use the names "ingoing" and "outgoing" to denote 
these paths below, even though in different coordinates they may not look ingoing 
or outgoing. 



We seek to promote our q = m test particles to black holes, expecting that they 
will follow "ingoing" or "outgoing" paths as in (4). Such black hole solutions should 
be quite complicated in static coordinates. For example, moving charged black 
holes will generate magnetic as well as electric fields. However, in cosmological 
coordinates the motion of such minimum energy particles is quite simple. 

The RNdS solutions can be rewritten in cosmological coordinates. For Q = M 
this has the form 



For M = this is just the standard cosmological form of the deSitter metric. 
We will call the coordinate system with H > expanding and that with H < 
contracting. The static Killing vector is given by 



The norm of the Killing vector vanishes at the horizons, which implies that the 
deSitter horizon r + and the outer black hole horizon r_ are located at Har±Vt 2 = 1, 



Cosmological Coordinates 




(1) 




(2) 
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or 

r± = 2^( 1 - 2M l ff l ± ^ 3 ^)' < 3 » 

There is also an inner black hole horizon at negative r in these coordinates. Note 
that the products a(t)r± are constants and also that the metric (1) is non-singular 
at r±. One can check that the surface gravity k at the two horizons, defined by 

JVa = "<a, (4) 

yields the correct value (4) for their common temperature. 

The transformation between the static and cosmological coordinates is given 

by 

a(t)r = R-M, t = T + h(R), h'{R) = - {R ^) V {R) (5) 

Integrating this transformation near the black hole and deSitter horizons, one finds 
that the expanding coordinates cover the past black hole horizon and the future 
deSitter horizon. Likewise, the contracting coordinates cover the past deSitter 
horizon and the future black hole horizon. Thus, the metric (1) with H > 
actually describes the white hole portion of the extended spacetime, while the 
metric with H < describe the black hole part. This can be confirmed by looking 
at null geodesies. The 2-sphere r = r_ is an outer trapped surface for H < 0. Null 
rays cannot get out. For H > it is an inner trapped surface. Null rays cannot 
get in. 

We can look at the paths of q = m test particles in (1). The paths (4) that were 
"outgoing" in static coordinates, stay at constant spatial comoving coordinate in 
the expanding coordinates, whereas they move out "doubly" fast in the contracting 
coordinates (and vice- versa for the "ingoing" test particles). 

It is useful in understanding the new solutions below to look at the motion of 
these particles in a little more detail. Consider an "outgoing" q = m test particle 
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in expanding coordinates (H > 0). It stays at a constant coordinate radius. The 
deSitter and white hole horizons r±, however, redshift like 1/a. At early times the 
horizons are both at large radii and the test particle is inside both. At some later 
time the white hole horizon sweeps past it, and then at some still later time the 
deSitter horizon sweeps past it. In static coordinates, we could follow the portion 
of the particles path between the two horizons. We would see it leave the white hole 
horizon at T = — oo and move out to the deSitter horizon at T = +00. Likewise, 
an "incoming" test particle, which stays at a fixed radial coordinate in contracting 
coordinates (H < 0), starts outside the past deSitter horizon, moves in through 
it and then through the black hole horizon. In static coordinates, this looks like 
a particle leaving the deSitter horizon at T = — 00 and getting to the black hole 
horizon at T = +00 



Cosmological MP Solutions 

The metric (1) for one \Q\ — M RNdS black hole in cosmological coordinates 
closely resembles the MP metric (1) for the case of one black hole. This suggests 
a simple form for the cosmological MP solutions 

ds 2 — —Q~ 2 dt 2 + a 2 (t)Q 2 (dx 2 + dy 2 + dz 2 ) , A t = n~ 1 , a(t) = e Ht , 
n=l + J2—, n = J{x - a*) 2 + (y- yi ) 2 + (z - *)2, H = ± 

(1) 

Indeed, we will show that more generally, a metric and gauge field of the form (1) 
solve the Einstein-Maxwell equations with cosmological constant A, if d t (a(t)Q) = 
a and V 2 fi = 0, where V 2 is the flat space Laplacian. 

First consider the constraint equations. The Hamiltonian, momentum and 
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Gauss' law constraints on a spatial slice are given by 

H = -^R + - (ir ij irij - ±tt 2 ) = -16np 

H k = - A(3) Vi7r ji = 16nJki (2) 

^WiE 1 = 0, 

where gij and 7C{j are the spatial metric and momentum, p and are the energy 
and current densities, and E l is the electric field. For a metric of the form (1), if 
we assume that dt(aVt) = a then the momentum is given by 

7Ty = -2-y/ggij. (3) 

CI 

Given this simple relation between 7i{j and the metric, the momentum constraint 
in (2) is satisfied with zero current. The 3-dimensional scalar curvature is 

where gij, Vj are the flat spatial metric and derivative. The energy density p has 
contributions from the Maxwell field and the cosmological fluid p = p ma x + Pcos- 
If A t = then 

8n Pmax = --L^Vi^Vj-a (5) 

From the expressions (3), (4) and (5), one sees that the Hamiltonian constraint can 
be satisfied by having the "black hole parts" and the "cosmological parts" vanish 
separately. That is, one has a solution if 

V 2 fi = 0, and Q =yp cos . (6) 

Further, the Gauss' law constraint is satisfied if (6) holds. Note that (6) can be 
satisfied by any time dependent, but spatially constant p cos . 
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Some insight into the nature of the solutions is gained from the matter equa- 
tions of motion, V a {T m b ax + T^ s ) = 0. The time component of this equation is 

Id a Id a , 

T^~nPcos "I - "-> — {Pcos ~t~ Pcos) 7^~T[Pmax 4" ^~Pmax = U (7) 

Hat a lildt a 

From (5) it follows that the Maxwell terms themselves sum to zero. Therefore, 
the part involving the cosmological fluid must also vanish independently. At this 
point there are two choices for a solution. From (6), p cos cannot depend on the 
spatial coordinates. If p cos is allowed to be time dependent, then, since Q has 
spatial dependence, the pressure p cos must be spatially dependent. This makes 
sense physically; ordinary matter would tend to accrete arround the black hole. 
The inhomogeneous pressure is needed to keep the density constant. In the one 
black hole case, such solutions have been studied by McVittie [11] and others 
(see e.g. [12]). The other choice, which we shall make, is to take p cos to be 
independent of time. We then have p cos = —p CO s = —A, the form for a cosmological 
constant. Again, this makes sense physically, because a cosmological constant 
cannot accrete*. Finally, we note that the full set of evolution equations for nij 
are straightforward to check and yield no more constraints. 



Geometric and Thermodynamic Properties 

The solutions (1) with H < appear to describe a system of "incoming" 
charged black holes. The solutions with H > would describe the time reversed 
situation; a system of "outgoing" charged white holes. The first thing to establish is 
that the solutions really do contain black hole (or white hole) horizons. One expects 
this to be the case. However, it is not straightforward to locate the horizons. The 
solutions, with more than one black hole, do not appear to have a stationary 
Killing vector. Thus, one cannot simply look for the surfaces on which the Killing 

* An idealized model of a cosmic string, spacetime minus a wedge, can be embedded in any 
flat Robertson- Walker background since the string does not accrete [13] 
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vector becomes null. One can look for apparent horizons (boundaries of regions of 
trapped surfaces) in a given spatial slice, but this too is complicated by the lack 
of symmetry. 

The situation does simplify for early and late times, when the holes are either 
far apart or close together. For concreteness, consider two "incoming" black holes 
in a background with H < 0. At early times the holes are far apart, and near 
each one the metric approaches that for a single hole. In this limit, one can verify 
that there are regions of outer trapped surfaces surrounding each of the holes. 
These regions extend out to radii r-j, given by (3). As time goes on, the universe 
contracts, and the coordinate size of each of these apparent horizons grows. Their 
shapes will distort due to the presence of the other hole. At some point in time 
the two apparent horizons will collide with one another and presumably merge. 
Indeed, in the late time limit one can verify that outer trapped surfaces surround 
both of the holes. 

To summarize, when H < 0, one can show that the two objects are first 
surrounded individually by outer trapped surface, and then later there is an outer 
trapped surface that surrounds both. For white holes with H > 0, the situation 
is time reversed. At early times inner trapped surfaces surround both objects 
together. Later on, the objects are surrounded only separately by regions of inner 
trapped surfaces. These results agree with what one expects from the area theorem. 
Black holes merge, and white holes split. 

Another question is whether an extension of the solutions with H > 0, which 
covers the black hole portions of the "outgoing" holes (and likewise for the white 
hole portions of the "ingoing" black holes with H < 0). For the case of one black 
hole, (5) gives an explicit coordinate transformation from expanding to contracting 
coordinates. It appears that making this transformation locally about one of the 
"outgoing" holes (as t — > oo), does indeed extend the spacetime to cover a black 
hole horizon. However, aspects of this transformation are still confusing; e.g. how 
the different regions fit together. We defer an explicit presentation to future work. 
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For the purposes of discussion in this paper, we will assume that the "outgoing" 
white hole spacetimes are extendable to coordinates that cover black hole horizons 
and past deSitter horizons. Likewise, the "incoming" black hole solutions are 
assumed to be fully extendable. A related question is whether the black holes must 
all be "outgoing" or all "incoming" , or whether there can be arbitrary combinations 
of "incoming" and "outgoing" holes. 

Although a full understanding of the thermodynamics of these solutions must 
await a better understanding of their analytic extensions and horizon structures, 
the thermodynamics appears to be quite interesting. Given that there is no station- 
ary Killing vector, the usual definition of temperature in terms of surface gravity 
at a horizon (2) does not work. Indeed, these solutions appear to be highly non- 
equilibrium. Still, as for the horizon structure, we can give a simple description of 
the thermodynamic behavior at early and late times, as follows. 

Recall that for one Q = M RNdS black hole, the deSitter and black hole hori- 
zons are in thermal equilibrium at a common temperature = T^s = T(M, A) 
given by (4). Now consider two black holes in a deSitter background, both sat- 
isfying Q = M, but with different masses Mi and M2. If the black holes arc 
widely enough separated to be out of causal contact, then each black hole will 
have its own distinct deSitter horizon. Also, in the region near each black hole 
there will be an approximate static Killing vector that can be used to define an 
approximate temperature. Each of the black holes will be in approximate thermal 
equilibrium with its own deSitter horizon at temperatures T(M\, A) and T(M2, A) 
respectively. If the two black holes are "ingoing" with respect to one another, then 
this will be the situation at very early times. Later the black holes come into causal 
contact and eventually merge into a single black hole with mass M = Mi + M2. 
At very late times, there will again be an approximate static Killing vector. The 
black hole horizon will be in equilibrium with the deSitter horizon at temperature 
T(Mi + M 2 ,A). 
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Conclusions and Further Questions 



We have presented solutions to the Einstein-Maxwell equations with a positive 
cosmological constant, which plausibly represent a collection of charged black holes 
moving either towards or away from one another. We have also given a description 
of the horizon structure and thermodynamics of these solutions in the early and 
late time limits. There is clearly a good deal more work to be done on these 
solutions; some of which has been noted in the previous section. 

Let's start with questions about the classical properties of the solutions. It is 
important to know whether these solutions can be fully analytically extended, as 
has been assumed above, and to have a clear picture of their horizon structures. 
It should be especially interesting to study the regime in which the black holes 
are merging (the white holes splitting). These solutions might give insight into 
an approximation which treats astrophysically interesting mergers of black hole 
binaries. Although, note that here the mergers take place without any gravitational 
radiation. 

The thermodynamics of the intermediate state, where the black holes are dis- 
tinct, but still in causal contact, should be interesting. Will the masses of the 
black holes change due to emission and absorption of Hawking radiation during 
this period? If so, then they should emerge from this nonequilibrium state with 
charges in general differing from their masses. Eventually then, each black hole 
should exchange radiation with the background until it again reaches a Q = M 
state. A first step towards understanding the exchange of energy, would be to 
study what a particle detector sees, if follows the path of a q = m test particle in 
a Q = M RNdS background. 

The splitting and merging of holes raises interesting questions about the pa- 
rameter space of A = solutions. Extreme RN black holes are regarded as solitons 
of general relativity, satisfying a kind of Bogomolnyi bound [14]. But in at least one 
respect they appear to differ from other solitons. Consider a magnetic monopole in 
the Bogomolnyi limit with two units of topological charge. The solution will have 
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zero modes corresponding to the possibility of breaking it up into two monopoles, 
each having unit topological charge. By analogy, we would expect an extreme RN 
black hole to have zero modes, corresponding to the possibility of breaking it up 
into smaller extreme black holes. By the area theorem, though, this cannot be the 
case. On the other hand, a single Q = M RNdS white hole appears to be unstable. 
It can be split into any number of charge equal mass fragments, which are then 
carried apart by the expansion of the universe. In the limit of zero cosmological 
constant, this may correspond to marginal stability of an extreme RN white hole. 
This picture is supported by analysis at the level of test particles. A q = m test 
particle can stay at rest in a Q = M RN background, whether it is inside or outside 
the event horizon. Hence, there should be analogues of the MP solutions describing 
merged black holes. 

Another interesting set of questions involves supersymmetry. Multi-object so- 
lutions are usually associated with Bogomolnyi bounds arising from an underlying 
supersymmetry of the solution. Romans [2] has noted, though, that the relevant 
supersymmetry (coming from N = 2 Yang-Mills supergravity) is consistent only 
with A < 0. Our solutions, then, are not supersymmetric, at least in this sense. 
On the other hand, for A < 0, while the Q = M RNdS holes are supersymmetric, 
they are also naked singularities. It should be interesting to see whether these can 
also be assembled into multi-hole solutions and to understand the role played by 
supersymmetry [15]. 



APPENDIX A 

In reference [9], Wald asked whether you can destroy a black hole by overcharg- 
ing it. One might think that, by throwing in particles with a high charge to mass 
ratio, one could charge the hole past the extremal limit of Q = M. It turns out [9] 
that for a charged particle to get over the Coulomb barrier into the hole, it must 
have more energy greater than or equal to its charge. 
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Here we do the analogous calculations for the extreme RNdS black holes [2]. 
For a metric of the form 



ds 2 = -V(r)dt 2 + V(r) _1 dr 2 + r 2 dQ 2 (1) 

with gauge field At(r) = —Q/r, the equation of motion for a test particle (of energy 
E, rest mass m and charge q) on a radial geodesic is given by (3) 

%) =-V(r) + -^(E + qA t (r)) 2 . (2) 

In order for the particle to get into the hole it must, at least, reach the event horizon. 
Looking for the minimum energy such particle, we set E = m and dr/dr = at 
the horizon. At the horizon radius p, we then have (V(p) = 0) 

i-^y>o. ( 3) 



This translates to 



mp 



™>Q-. (4) 
q p 

For the extreme RNdS black holes we have [2], 

M = p(l-hp 2 ), Q 2 = p 2 (l-Ap 2 ). (5) 

From these we can compute how the mass and charge of an extremal hole change 
with each other. We find 

dM r Q 

which coincides with the bound (4). Hence an extreme RNdS black hole cannot 
be pushed over the limit. Given the nontrivial relation (5) between the charge and 
mass of the extreme black holes, this is a somewhat striking confirmation of cosmic 
censorship. 
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